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J /if) production in p+p collisions at y/s = 200 GeV has been measured in the PHENIX ex- 
periment at the Relativistic Heavy Ion Collider (RHIC) over a rapidity range of —2.2 < y < 2.2 
and a transverse momentum range of < pr < 9 GeV/c. The statistics available allow a de- 
tailed measurement of both the pr and rapidity distributions and are sufficient to constrain pro- 
duction models. The total cross section times branching ratio determined for J/ip production is 
Bu ■ o 3 J* = 178±3 stat ±53 sys ±18 norm nb. 

PACS numbers: PACS numbers: 



J/ip are produced in hadronic collisions involving hard 
processes that proceed primarily through diagrams in- 
volving gluons, such as gluon-gluon fusion. Once the cc 
pair is produced it must evolve through a hadronization 
process to form a physical J/ip. While this production 
has been extensively studied, the details of the produc- 
tion mechanism and hadronization remain an open ques- 
tion. Attempts at a consistent theoretical description 
of J/ip production have been made, but it has proven 
difficult to reproduce both the observed cross sections 
and polarization [l|, 0, 0, 0] • An additional complication 
is that nearly 30-40% of the measured J/ip yield results 
from feeddown of higher mass states (ip\ Xc) M, reducing 
the observed polarization with respect to that expected 
from directly produced J/ip. 

The color-singlet model [6J, which generates a color 
singlet cc pair in the same quantum state as the J/ip, 
underpredicts the measured J/ip cross section by approx- 
imately an order of magnitude Q. Alternatively, the 
color-octet model includes color octet cc pairs that radi- 
ate soft gluons during J/ip formation Q. However, the 
color octet matrix elements are not universal Q and the 
predicted transverse J/ip polarization at high pt is not 
seen in the data 0, [j]. The color evaporation model, 
a more phenomenological approach, forms the different 
charmonium states in proportions determined from ex- 
perimental data for any cc pair that has a mass below 
the DD threshold and predicts no polarization. Finally, 
a recent perturbative QCD calculation including 3-gluon 
diagrams is able to successfully reproduce both the ob- 
served cross section and polarization results 

A fundamental understanding of the J / ip production 
process is also critical to defining the configuration of the 
produced cc state since this will have direct implications 
on the interaction of this state with both cold nuclear 
matter in proton or deuteron-nucleus collisions and with 
the high-density partonic matter observed in high-energy 
heavy-ion collisions. High quality experimental results 
over wide kinematical and collision energy ranges are re- 
quired to constrain models and to provide an improved 
understanding of J/ip (and other heavy quarkonia) pro- 
duction. 

In this paper, J/ip production in p+p collisions at 
y/s = 200 GeV measured by the PHENIX experiment 
at RHIC is reported. The J/ip cross section and trans- 
verse momentum distributions arc studied in the mid 
(\y\ < 0.35) and forward (1.2 < |y| < 2.2) rapidity re- 



gions. The data presented were collected during the 2005 
RHIC run and exceed by more than one order of magni- 
tude the previously reported number of J/ip [3, 11 1 . 

A detailed description of the PHENIX experiment is 
provided in [T3 |. At mid-rapidity the drift chambers 
(DC) , ring imaging Cerenkov detectors (RICH) , and elec- 
tromagnetic calorimeters (EMCal) are used to detect 
J/ip — > e + e~ decays in two arms each covering A(p = 90° 
in azimuth. The muon detectors, consisting of cath- 
ode strip tracking chambers in a magnetic field (MuTr) 
and alternating layers of steel absorber and Iarocci tube 
planes (MuID), are used to measure J/ip — * /i + /J~ at 
forward and backward rapidities over A<p — 360°. 

The data were recorded using a minimum bias trig- 
ger that requires at least one hit in each of the beam- 
beam counters (BBC) at forward and backward rapidity, 
3.0 < |tj| < 3.9. Di-electron events must pass an addi- 
tional trigger that consisted of an OR between the level- 1 
electron and photon triggers. The electron trigger re- 
quires matching hits between the EMCal and RICH in a 
small angular area with a minimum energy deposition of 
0.4 GeV in any 2x2 patch of EMCal towers. The photon 
trigger requires a minimum energy deposition of 1.4 GeV 
in any 4 x 4 set of overlapping EMCal towers. A J/ip trig- 
ger efficiency of 96% was achieved within the vertex range 
\%vtx\ < 30 cm. Di-muon triggered events were selected 
using an online level- 1 trigger that requires at least two 
particles penetrate the MuID. One particle must pene- 
trate the entire MuID while the second has a minimum 
penetration depth of 3 out of the 5 pairs of detector and 
absorber planes. Approximately 92% of the J ftp's within 
l^tel < 30 cm fulfill this requirement. As part of the re- 
construction chain a level-2 filter, which consists of a fast 
reconstruction of the particle trajectory in the MuTr and 
MuID is applied. Events are accepted by this filter when 
at least two particles penetrate the entire MuID and have 
a reconstructed invariant mass > 2.0 GeV/c 2 . After ap- 
plying cuts on the collision vertex position and quality 
assurance criteria, the sampled statistics corresponds to 
2.6 pb^ 1 in the dielectron analysis, 2.7 pb" 1 in the muon 
arm covering 1.2 < y < 2.2 and 3.5 pb _1 in the muon 
arm covering —2.2 < y < —1.2. 

At mid-rapidity electron candidates are charged tracks 
associated with at least two hit phototubes in the RICH 
and one EMCal hit with a position matching of ±4 stan- 
dard deviations (er). The energy-momentum matching 
requirement is (E/p— 1)/ >-4.0cr. The number of J/ip 
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FIG. 1: Invariant mass spectra for a)J/tp — » e + e at |y| < 
0.35 and b) J/tp -> fi + fi~ at 1.2 < \y\ < 2.2 with the func- 
tional forms used to extract the number of reconstructed J/tp. 



is obtained by counting the unlike-sign dielectron pairs in 
a fixed mass window after subtracting the like-sign pairs. 
Figure [TJa) shows the invariant mass spectrum for di- 
electron pairs after subtracting the like-sign background. 
The mass window for counting the J/tp signal is 2.7 - 
3.5 GeV/c 2 or 2.6 - 3.6 GeV/c 2 depending on the number 
of the DC hits used to reconstruct the track. The J/tp 
counts are corrected for the continuum yield, which orig- 
inates primarily from open charm pairs and Drell-Yan, 
inside the mass window (10% ± 5%) and the fraction of 
J/tp outside of the mass window (7.2%±1.0%). Approxi- 
mately 1500 J/tp — ► e + e~ are obtained. The solid line in 
the figure is the sum of the J / tp line shape (dashed curve) 
and an exponential function (dot-dashed curve) describ- 
ing the continuum component. The J/tp line shape func- 
tion includes the detector resolution, the internal radia- 



tive effect [13| , and the external radiative effect evaluated 
using a GEANT [3] simulation of PHENIX detector. 

Muon track candidates are selected based upon their 
penetration depth in the MuID and the reconstructed 
track quality within the MuID and MuTr. The particle 
trajectory must contain at least 8 of 10 possible hits in 
the MuID and the position matching between the MuID 
and MuTr must be within 15(20) cm at positive (nega- 
tive) rapidity. The J/tp yield is obtained from the unlike 
sign dimuon invariant mass distribution by subtracting 
the combinatorial background estimated using an event 
mixing technique. Three functions, shown in Fig. [ljb), 
arc used to extract the J/tp yield. Single Gaussian + ex- 
ponential and double Gaussian + exponential functions 
are used to fit the J/tp peak, while the contribution from 
the physical continuum and background is estimated us- 
ing an exponential fit. The reported number of J/tp rep- 
resents the average of the fit values. A total of 8000 
J/tp — ^ are obtained. The systematic error asso- 

ciated with the signal extraction is estimated from the 



variation between the fits. 

The J/tp cross section in a given rapidity and trans- 
verse momentum bin is calculated according to, 
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2ttpt dydpT 2n p T Ap T Ay LA e roc e tr i g e^c ' 

where Bu is the J/tp dilepton branching ratio, Nj/^ is 
the measured J/tp yield, C is the integrated luminosity 
recorded by the minimum bias trigger, Ae roc represents 
the geometrical acceptance and reconstruction efficiency, 
and e t rig is the trigger efficiency. is the minimum 

bias trigger efficiency for events containing a J/tp and 
was determined to be 0.79 ± 0.02 [lfj. The cross section 
sampled by the BBC trigger, af£ x e^| c = 23.0±2.2mb, 
was used to determine the integrated luminosity. 

The Ae ICC and etrig terms are determined individually 
for the central arm and each muon arm based upon the 
detection of simulated J / tp processed using the real data 
analysis chain. Decay events are generated and propa- 
gated through a full GEANT simulation of the detector, 
which includes the specific details of the detector perfor- 
mance including the MuTr and MuID alignment, disabled 
anodes and MuID efficiency. For the dielectron analy- 
sis, corrections to account for the detector dead channel 
map, energy calibration and run-to-run variations in the 
detector active area were determined from single electron 
yields. The J/tp trigger efficiency is incorporated via a 
level- 1 trigger emulator tuned to describe the experimen- 
tal trigger response. For the dimuon analysis, the level-2 
filtering algorithms are applied to the simulated events. 
After reconstruction, the number of detected J / tp is com- 
pared to the number of simulated J/tp in a given rapidity 
and transverse momentum bin to determine the appro- 
priate correction factors. 

The systematic error associated with the measurement 
of the J/tp cross section can be divided into three cate- 
gories based upon the effect each error has on the mea- 
sured results. All errors are reported as standard de- 
viations. Point-to-point uncorrelated errors, such as the 
signal extraction systematic, which is bin dependent with 
typical values of 4% (5%) in the dimuon (dielectron) data, 
allow the data points to move independently with respect 
to one another. Point-to-point correlated errors allow the 
data points to move coherently within the quoted value. 
Their values amount to 10% (8%) for the detector ac- 
ceptance, 8% (4%) for the run-to-run variation in the 
detector efficiency, 4% (2.5%) for the J/tp transverse mo- 
mentum and vertex distributions, and 2% (2%) for the 
hardware efficiency of the detector. Finally, global sys- 
tematic errors allow the data points to move together. 
The dominant source of this error originates from the es- 
timation of the BBC triggering efficiency for minimum 
bias events, 9.7%, with additional contribution from the 
uncertainty in the estimation of the number of sampled 
minimum bias events, 1%, and e^A^' 2.5%. 
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Figure[2la) shows the transverse momentum spectra at 
both mid and forward rapidities, which are fit with the 
function, Ax (1 + (p T /B) 2 )" 6 to extract the (p|>. 
At mid-rapidity (p|> = 4.14 ± 0.18±g;|g (GeV/c) 2 and 
the x 2 P er degree of freedom (x 2 /ndf) is 23/19. At for- 
ward rapidity (p\) = 3.59 ± 0.06 ± 0.16 (GeV/c) 2 and 
the % 2 /ndf is 28/17. The first error is statistical and the 
second includes the systematic error from the maximum 
shape deviation permitted by the point-to-point corre- 
lated errors and from allowing the exponent of the fit 
function to be a free parameter. Although good agree- 
ment is found with the rapidity distribution and total 
cross section, previously published results [ill ] yielded a 
significantly lower (p T ) at forward rapidity than found 
here, even accounting for the quoted statistical and sys- 
tematic uncertainties. The increased statistics of the 
present data set allow for an improved understanding of 
the shape of the spectrum at forward rapidity due to 
the extended range in px and the finer binning at low 
Pt- The previous results have been revisited and it was 
found that the systematic error was underestimated. 

Figure [Hh) shows the ratio of the invariant cross sec- 
tion at forward and mid-rapidity. The ratio falls with pt 
and reaches a minimum of 0.5 above a px of 2 GeV/c. 
The data indicates that the forward rapidity pt distri- 
bution is substantially softer than mid-rapidity. This is 
attributed to the increase in the longitudinal momentum 
at forward rapidity leaving less energy available in the 
transverse direction. 



Global scale uncertainty: 10.1% 
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FIG. 2: (a) The mid and forward rapidity J/ip differential 
cross section times dilepton branching ratio versus pr and (b) 
the ratio of the mid and forward pr spectra. The vertical error 
bars are the statistical and point-to-point uncorrelated error 
and the boxes are the point-to-point correlated systematic 
error. The solid lines are the fits described in the text. 

The observed pt distributions are substantially harder 



than those for lower energy p+p and p+A collisions as 
expected from the increased phase space at higher en- 
ergy. Figure [3] shows the energy dependence of the av- 
erage (p T ) including lower energy points from the Super 
Proton Synchrotron, and Fermilab fixed target and Teva- 
tron measurements. A linear fit versus the log of the cen- 
ter of mass energy describes the general trends, although 
some variation is expected due to the differing rapidity 
ranges of the measurements and the use of p + A data for 
some of the points. 
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FIG. 3: PHENIX (p^) measurements compared to measure- 
ments at other energies 0, ["Hil . as a function of collision 
energy for J/ip production in p+p or p+A collisions. Also 
shown is a linear fit vs \n(y/s). 
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FIG. 4: The J/tp differential cross section times dilepton 
branching ratio plotted versus rapidity. The vertical error 
bars represent the statistical and point-to-point uncorrelated 
error and the boxes represent the point-to-point correlated 
systematic error. The curves are described in the text. 



Figured] shows the J/tp differential cross section vs ra- 
pidity. The statistics of these results are large enough 
to allow eleven rapidity bins to be measured compared 
to the five bins used in the previous measurement [ll|. 
Also shown are several model fits to the data. The dashed 



curve is a non-relativistic QCD calculation 18]. The dot- 
dash curve is a pQCD calculation that includes diagrams 
describing a third gluon, which is necessary to neutral- 
ize the otherwise colored states 0. This model fails to 
reproduce the steeply falling cross section observed in 
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the present data at large rapidity. An empirical dou- 
ble Gaussian fit (dot-dot curve) is able to reproduce the 
data best, but has no strong theoretical foundation. The 
data slightly favor a flatter distribution over the rapidity 
range \y\ < 1.5 than most models, but since the system- 
atic error on the mid and forward rapidity points are 
independent, a narrower distribution is not excluded. 

To determine the total cross section, the rapidity dis- 
tribution was fit with many theoretical and phenomo- 
logical shapes, including those shown in Fig. |U We 
obtain a total cross section times branching ratio of 
B u ■ a J v ^ = 178±3 stat ±53 sys ±18 norm nb. The absolute 
normalization error (norm) represents the uncertainty in 
the BBC trigger cross section. The systematic error (sys) 
is estimated from the maximum variation allowed by 
shifting the mid and forward rapidity data independently 
by their point-to-point correlated systematic errors. This 
result is consistent with our previous measurement [llj ]. 

We have presented J/ip results for p+p collisions at ^fs 
— 200 GeV that extend the reach in transverse momen- 
tum to 9 GeV I c. The measured px spectrum is harder 
than that observed at lower energies and also shows a 
softening at forward rapidity. The rapidity shape falls 
steeply at forward rapidity and can not be reproduced 
by the pQCD calculation in Q. Futhermore, the data 
slightly favor a flatter rapidity distribution than most 
models, but a narrower distribution is not excluded. 
These data not only constrain production models for 
heavy quarkonia, but also provide a critical baseline for 
similar studies in deuteron-nucleus and heavy-ion colli- 



OTKA-MTA, and US-Israel BSF. 



We thank the staff of the Collider-Accelerator and 
Physics Departments at BNL for their vital contribu- 
tions. We acknowledge support from the Department of 
Energy and NSF (U.S.A.), MEXT and JSPS (Japan), 
CNPq and FAPESP (Brazil), NSFC (China), MSMT 
(Czech Republic), IN2P3/CNRS, and CEA (France), 
BMBF, DA AD, and AvH (Germany), OTKA (Hungary), 
DAE (India), ISF (Israel), KRF and KOSEF (Korea), 
MES, RAS, and FAAE (Russia), VR and KAW (Swe- 
den), U.S. CRDF for the FSU, US-Hungarian NSF- 



t 
[1] 

[2] 

[3] 
[4] 
[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 
[13] 
[14] 

[15] 

[16] 

[17] 

[18] 
[19] 
[20] 



Deceased 

PHENIX Spokesperson: zajc@nevis.columbia.edu 

M. Beneke, M. Kraemer, Phys. Rev. D55, 5269 (1997); 

M. Beneke, I. Z. Rothstein, Phys. Rev. D54, 2005 (1996). 

F. Abe et al, Phys. Rev. Lett. 69, 3704 (1992); F. Abe 
et al, Phys. Rev. D66, 092001 (2002); S. Abachi et al, 
Phys. Lett. B370, 239 (1996); B. Abbot et al, Phys. 
Rev. Lett. 82, 35 (1999). 

R. Vogt, Phys. Repts. 310, 197 (1999). 

M. H. Schub et al, Phys. Rev. D52, 1307 (1995). 

I. Abt et al, Phys. Lett. B561, 61 (2003). 

E. L. Berger et al, Phys. Rev. D23, 1521 (1981); R. Baier 

et al, Phys. Lett. B102, 364 (1981). 

G. T. Bodwin et al, Phys. Rev. D51, 1125 (1995); erra- 
tum Phys. Rev. D55, 5853 (1997). 

J. K. Mizukoshi, SLAC-PUB-8296, |hep-ph/991i384| 
(1999). 

V. A. Khoze et al, Eur. Phys. J. C39, 163-171 (2005) 
and private communications. 

S. S. Adler et al, Phys. Rev. Lett. 92, 051802 (2004). 
S. S. Adler et al, Phys. Rev. Lett. 96, 012304 (2006). 
K. Adcox et al , Nucl. Instr. Me th. A499, 469 (2003). 
A. Spiridonov, |hep-ex/0510076| (2004) . 
GEANT 3.2.1, CERN Computing Library, 



http://wwwasd.web.cern.ch/wwwasd/geant/index.html 



J. K. Yoh et al, Phys. Rev. Lett. 41, 684 (1978); and 
private communication.. 

O. Drapier, These d'habilitation, Universite 

Claude Bernard - Lyon 1, 1998; available at 

http://na50.web.cern.ch/NA50/theses.html 

D. Acosta et al, Phys. Rev. D71, 032001 (2005). The 

(px) was extracted using a fit of the form A x (1 + 

(px/B) 2 )- 6 . 

F. Cooper et al, Phys. Rev. Lett. 93, 171801 (2004); and 
private communication. 

S. S. Adler et al, (Au+Au J/i/j paper, preprint number 
soon). 

S. S. Adler et al, (Cu+Cu J/if) paper, preprint number 
soon). 



